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Abstract 

We construct off-diagonal parton distributions defined on the interval < X < 1 
starting from the off-forward distributions defined by Ji. We emphasize the partic- 
ular role played by the symmetry relations in the "ERBL-like" region. We find the 
evolution equations for the off-diagonal distributions which conserve these symme- 
tries. We present numerical results of the evolution, and verify that the analytic 
asymptotic forms of the parton distributions are reproduced. We also compare the 
constructed off-diagonal distributions with the non-forward distributions defined 
by Radyushkin and comment on the singularity structure of the basic amplitude 
written in terms of the off-diagonal distributions. 



1 On leave from H. Niewodniczanski Institute of Nuclear Physics, ul. Radzikowskiego 152, Krakow, 
Poland. 



1 Introduction 



It is well known that the cross section of hard scattering processes (such as deep inelastic 
scattering, the production of large pt jets, etc.) can be written as the sum of parton distri- 
butions multiplied by the cross sections of hard subprocesses calculated at the parton level 
using perturbative QCD. That is we can factor off the long distance (non-perturbative) 
effects into universal, process independent, parton distributions (fi(X, fi 2 ) with i = q,q, g) 
specific to the incoming hadrons. X is the longitudinal fraction of the hadron's momen- 
tum that is carried by the parton and fi is a scale typical of the hard subprocess. The 
parton distributions are given by the matrix elements (P\0\P) where O is a twist-2 quark 
or gluon operator, and P represents the full set of quantum numbers of the hadron. To 
be specific we will be concerned with a proton taking part in unpolarised reactions. Thus 
P will represent the 4-momentum of the proton. 

Calculating the parton distributions from first principles is one of the most challenging 
problems in non-perturbative QCD. The most promising approach is lattice QCD, but 
much remains to be done. On the other hand, from a practical viewpoint, the parton 
distributions of the proton are determined with good precision from global analyses of deep 
inelastic and related hard scattering data. The distributions fi{X,jj?) are parametrized 
as a function of X at some starting scale //q and then evolved using the DGLAP equations 
of perturbative QCD to higher /i 2 values relevant to the data to be fitted. 

Recently there has been much interest in off-diagonal (also called off-forward 

by Ji [JTJ] or non- forward by Radyushkin |3J) distributions which are given by matrix ele- 
ments (P'\0\P) in which the momentum of the outgoing proton is not the same as that 
of the incoming proton. For example, the amplitudes for processes such as deeply virtual 
Compton scattering (7*^ — > jp) or vector particle electroproduction (7*^ — > Zp or J/ipp) 
depend on off-diagonal distributions. Since P 7^ P' the parton returning to the proton 
has a different momentum to the one which is outgoing, and so we need two momentum 
variables to specify the off-diagonal distributions. The Ji and Radyushkin distributions, 
which are denoted by if (2, £) and J-^{X) respectively, differ in their choice of the defining 
four vector. Ji chooses the momentum fractions x and £ with respect to the average of the 
incoming and outgoing proton momenta P — \{P + P'), whereas Radyushkin defines X 
and ( with respect to the incoming proton momentum P. The former has the important 
advantage that it is easier to impose the symmetry requirements, while the latter has the 
advantage that it is close to the definition used for the conventional (diagonal) distribu- 
tions. Our aim is to clarify the relation between the two formulations. We find that they 
are not equivalent unless specific conditions are imposed on Radyushkin's non-forward 
distributions. We show this by a direct construction of distributions defined in the range 
< X < 1 which are equivalent to Ji's off-forward distributions. 

Let us neglect, for the moment, the gluon distribution. The quark distribution 
H q (x,£,), defined by Ji, covers the interval — 1 < x < 1 and generates two distinct distribu- 
tions which we denote?] by T q {X, () and Fq(X, () with < X < 1. Over the region X > ( 

2 For the reasons given below we must use a notation which distinguishes between the distributions 
F(X, C) constructed from H and the non-forward distributions JF^(X) defined by Radyushkin. 
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the two functions T q and Tq are independent. On the other hand in the region X < ( they 
are related to each other, with the consequence that the non-singlet and singlet combina- 
tions possess a symmetry about X = (/2. We obtain evolution equations for jF starting 
from the evolution equations for the off-forward distributions H. We find that they differ 
from the evolution equations for the non-forward distributions |5], [| by additional terms 
which are essential to preserve the symmetry properties in the ERBL-like region. We 
also found that the basic amplitude for deeply virtual Compton scattering (DVCS) has a 
different singularity structure to that given by the non-forward distributions T . 

The outline of the paper is as follows. To establish notation we quickly review in 
Section 2 the salient features of the conventional (diagonal) parton distributions H(x) 
with support — 1 < x < 1. Section 3 reviews the extension of these ideas to the off- 
diagonal distributions if(x,£) that were introduced by Ji [|TJ. In Section 4 we transform 
the distributions if(x,£) into distributions TiX, Q with < X < 1, and demonstrate 
that T must satisfy symmetry relations for X < In Section 5 we give the evolution 
equations for the ^(X, () and present numerical solutions. The complete form of the 
evolution equations is given in the Appendix. In Section 6 we discuss the relation between 
the distributions T and the non-forward distributions T of Radyushkin. In the same spirit 
we discuss the differences in the singularity structure of the DVCS amplitude. Finally 
Section 7 contains our conclusions. 



2 Conventional parton distributions 

In order to introduce off-diagonal distributions it is most convenient to first recall the 
definition of the conventional (diagonal) parton distributions in terms of light-cone coor- 
dinates {x ± = (x° x 1 , x 2 ) and in the light-cone gauge (A + = 0) JjlJ. For instance 
the quark distribution H q (x) is given in terms of the matrix element of a light-cone bilocal 
operator 

Note that the matrix element is diagonal in the four momentum P of the proton. For 
simplicity we do not show either here, or throughout the paper, the renormalization scale 
dependence of H q and of the other parton distributions that we discuss. 

To see the parton content of the distribution H q we make a Fourier expansion of 
(the light-cone-plus or 'good' component) ip + of the quark field, in terms of the quark 
annihilation operator b and the antiquark creation operator S. Similarly ip + is expanded 
in terms of and d, and then the integration over y~ in ([l]) is performed. It is found 
that H q is only non-vanishing in the interval — 1 < x < 1 with the term tfb contributing 
for x > and contributing for x < |TT| 

H q {x) = ^ J ^^^[(P\b{{xP\k T )bixP\k T )\P)9{x) 

(2) 

- (P\d{(-xP + , k T ) d x {-xP + , k T )\P) 6{-x) 
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where A is the helicity of the quarks. The tfb term corresponds to the emission of a 
quark (carrying a fraction x of the proton's momentum) and its subsequent reabsorption 
within the proton. Similarly the cftd contribution describes the emission and subsequent 
reabsorption of an antiquark. The two possibilities are sketched in Fig. 1. Thus the single 
distribution H q with support in the interval — 1 < x < 1 embodies both the familiar q 
and q distributions, defined on the interval < x < 1, which thus are identified with the 
two terms accompanying the theta functions in (Q) in the following way 

{q(x) for x > 
(3) 
—q(—x) for x < 0. 

We may form the valence and singlet quark distributions in terms of H q 

q(x) - q(x) = H q (x) + H q (-x) = Hj(x) 



(4) 



E + <?»] = E [H q {x) - H q (-x)} = H s (x) 



where the sum is over the quark flavours. Clearly over the full interval — 1 < x < 1 the 
valence and singlet quark distributions satisfy the symmetry relations 

H V q {x) = H v q {-x) 

(5) 

H s (x) = -H s (-x) . 

In a similar way we may introduce H g (x) = xg(x) where g(x) is the familiar gluon 
distribution. In the light-cone gauge 

H 9 (x) = ± j ^e-™ p+ y- (P\F+»(0,y-,0)F+(0,0 7 0)\P) , (6) 

where F^ v is the gluon field strength tensor and where the summation over the colour 
label has been suppressed. Due to Bose symmetry we have 

H a (x) = H g {-x) . (7) 



3 Off-diagonal distributions 



The distributions H q introduced in ([]]) may be generalized to allow for matrix elements 
which are off-diagonal in the four momentum of the proton 

H q {x,£,t) = \J ^-e- p+ y- (P'\M0,y-/2,O) I 7 +^(0, -y"/2, 0)|P) , (8) 

where we consider only the distributions which conserve the proton helicity and which 
describe unpolarized quarks. Since A = P — P' ^ the distribution H q (x,l;,t) now 
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contains two extra scalar variables, in addition to the Bjorken x variable. The variable t 
is the usual t-channel invariant, t = A 2 , and the variable £ is defined by 

6P + = ~A + , (9) 

where P = \{P + P'). This choice of variables^ is due to Ji 0-0 and enables symmetry 
to be imposed between the incoming and outgoing proton. That is Ji uses the symmetric 
combination P of their momenta as the defining direction, and calls the H q off-forward 
distributions. The distributions H q are real, and the symmetric choice of variables has 
the considerable advantage that, due to time-reversal invariance and hermiticity, the dis- 
tributions are even functions of £ 0] 

H q (x,£,t) = H q {x,-£,t) . (10) 

Since we will perform our analysis for fixed t, concentrating on the x and £ dependence, 
we shall omit the t dependence from now on. 

To see the physical content of the off-diagonal distributions H q we again Fourier expand 
ip and ip in terms of the quark creation and annihilation operators. Since the distributions 
are even in £ we may take £ > 0. In this way we obtain the generalization of eq. (Q) || 

TT , . 1 f d 2 k T <N 

HJx.g) = / — . (11) 

qV > q) 2P+J 2^^(271)3 1 ' 

J2 [{P'\b{ ((x - £)P+, k T - A T ) b x ((x + 0P + , k T )\P) 9(x > 

A 

+ (P'\d x ((-x + C)P + , -k T + A T ) 6_ a ((x + £)P + , k T )\P) 6{-Z <x<0 

- (P'\4 {{-x - 0P + , k T - A T ) d x ((-x + 0P + , k T )\P) 9{x < -£)" . 

Fig. 2 gives a pictorial description of the content of dTlf) . Diagrams (a) and (c), which 
arise from the tfb and cftd terms in x/ji/j, generalize Figs. 1(a) and (b) respectively. For 
example the first diagram corresponds to the emission of a quark of momentum k from 
the proton followed by its absorption with momentum k — A. Thus for x > £ and x < — £ 
the off-diagonal distribution H q generalizes the familiar quark and antiquark distributions 
and will evolve according to modified DGLAP equations. Diagram (b), corresponding to 
the middle region, — £ < x < £, does not have a counterpart in Fig. 1. This diagram, 
which arises from the db term in ipip, corresponds to the emission of a quark-antiquark 
pair. In this region H q is a generalization of the proton form factor and will evolve 
according to modified ERBL equations Thus in this domain H q may be regarded as 



a generalization of the probability distribution amplitude which occurs in hard exclusive 
processes. 

Just as for the diagonal case, we introduce valence and singlet quark distributions 
analogous to (f|) 

<(x,£) = H q {x,£) + H q (-x,0 = <(-x,£) , (12) 



3 Note that Ji defines A = P' - P. 
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H S (x, = E [H q (x, - 0] = -^ 5 (-a:, • (13) 

Thus in addition to the symmetry under £ — > — £, the distributions have symmetry or 
antisymmetry under x — > — x. Also, in analogy to (0), the off-diagonal gluon distribution 
satisfies 

ff p (x,0 = H g (-x,£). (14) 
are identical to those introduced by Ji [§, §]0 except that 

Pf s (x,0 = xHf(x,0- (15) 

On account of the extra factor x, the gluon distribution (|T5|) is not required to be zero at 
x = 0, unlike the situation for H^ 1 (see also for a relevant discussion). 



The distributions (Pj)-(i4 



4 Off-diagonal distributions on the interval [0, 1] 

So far we have considered the off-diagonal distributions H q (x,£), introduced by Ji |], ||, 
and defined on the interval — 1 < x < 1. As noted above P + P' is taken as the defining 
direction, so that symmetry is imposed between the incoming (P) and outgoing (P') 
proton momenta. This variable £ was defined in (|]) by 

A = (P — P') — £(P + P') , (16) 

where for simplicity we have omitted the light-cone plus superscript (see @). 

To make direct contact with conventional partons we may introduce alternative off- 
diagonal distributions J- q (X, £) defined on the interval < X < 1 such that the initial 
parton carries a positive fraction X of the proton's longitudinal momentum. That is we 
take P as the defining direction. Thus the counterpart to ( [IB] ) is 

A = CP (17) 

with < £ < 1. This is exactly analogous to the approach introduced by Radyushkin f|, |j 
in the construction of the non- forward distributions J-q{X). However our construction of 
the distributions F q (X, £) presented below is different to that of M H]. From (pip and 
flT?D it follows that 

f = """T • (18) 

4.1 The relation between the distributions H and T 

In this subsection we first define the off-diagonal distributions J- q {X, () with X in the 
interval [0,1] starting from Ji's distributions H q (x,£) with x in the range [-1,1]. Then we 
explore the symmetry relations satisfied by the fi q {X, (). 

4 Note that in going from Ref. to Ref. |^| Ji has redefined £/2 by ^. 
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If we compare the momentum fraction carried by the emitted parton in Fig. 3 with 
those in Figs. 2(a) and 2(c), then we see that two different transformations are relevant in 
reducing the interval — 1 < x < 1 covered by H q {x,£) to the interval < X < 1 covered 
by J- q (X, £). First, from Fig. 2(a), we have the transformation 

• Y ' = TTT • < 19) 

which takes the interval x% G [— £, 1] into X\ G [0,1]. Simultaneously £ is transformed 
into (. Secondly, from Fig. 2(c), we have the transformation 

*» = l=f . < 2 °) 

which takes £2 G [ — 1,6] into X2 G [0, 1]. Now, — £ is transformed into (. In this way we 
introduce two distinct off-diagonal distributions jF q and 

(21) 

where £ = C/(2 — C) an d the inverse relations 

*i ~ C/2 C/2 ~ X 2 

Xl = T3c72" ' X2 = T3c72" (22) 

follow from (|18| -^0|). We stress that as X 12 cover the range [0, 1], the corresponding X\ 
and X2 cover respectively the ranges [— £, 1] and [— 1,£], as shown schematically in Fig. 4. 
The factors ±(1 — £/2) _1 in (|2lf ) arise from the translation of the measure fix to dX. 

In the limit that ( (and £) — > we have from (|3|) 

^ g (X,0) = # 9 (X,0) = g(X) 

(23) 

J^(X,0) = -^(-X,0) = g(X), 

which is an additional motivation for using the quark and antiquark subscripts to differ- 
entiate between the two functions jF q and jF q . 

Finally, due to the symmetry relation ([H]), the gluon distribution may be defined in 
the range < X < 1 by either of the transformations (p2|). That is we have 



r,(X,Q = ^H.fc&.i) = jJ p ir.(^=|,£l. (24) 
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4.2 Symmetry relations 



From Fig. 4 we see that in the DGLAP-type regions (x > £ or x < — £) H q is transformed 
respectively into independent functions Fq(X) and Fq(X) with X > £. On the other hand 
in the ERBL-type region (— £ < x < £) the distribution if 9 generates functions •Fg(-X') 
and Pq{X) with X < £ which are no longer independent. Indeed for X < £ we have 

£^-;n 1 frf C-*-C/2 \ 1 h( ^ 2 ' x \ 

^ Aj i - C/2 9 V i-C/2 J " 1-C/2 H1-C/2J 

= -^-W, (25) 

where for simplicity we do not indicate the additional explicit £ or £ dependence of the 
distributions. 

Eq. fl25|) is the basic symmetry relation for the off-diagonal quark distributions which 
indicates that in the ERBL-like region the quark and antiquark distributions are not 
independent, unlike the case in the DGLAP-like region. The physical reason for this can 
easily be understood by looking at Fig. 2b. In the ERBL-like region we can define the 
off-diagonal distributions with respect to the first emitted parton being either the quark 
with momentum x + £ or the antiquark with momentum £ — x. The latter possibility 
corresponds to the exchange of the annihilation operators in eq. (pTTj) , which is the origin 
of the — sign in relation (E5p. 

We may form the non-singlet and singlet combinations of the quark and antiquark 
distributions. From (O) we have 



(26) 



&{X) = -——HY[—%-\ = F q (X) - ^(X) , 




which in the region X < ( satisfy symmetry relations resulting from ([25] 

^(C-x) = ^(x), 

f s ((-X) = -F S {X). 



(27) 



It is straightforward to show for X < Q that the gluon distribution (^J) satisfies a similar 
relation 

f g {C-X) = T g {X). (28) 

These properties are well-illustrated by Fig. 5. The upper plot shows an example of the 
off-diagonal distribution H q (x, £) with x G [—1, 1], for £ = 0.5. The middle plot shows the 
transformation of this distribution into the two functions F q (X, () and T q {X, () of (|2lf ) 
with X G [0, 1]. Their behaviour shows that the symmetry relation (|25| ) is clearly satisfied 
in the region < X < £. Finally, the lower plot shows the behaviour of the non-singlet 
TY and the singlet T s combinations. The symmetry of TY and antisymmetry of JF 5 , 
about the point X = (/2, are clearly evident in the region < X < £. 
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5 Evolution equations 



Just as we constructed jF directly from the off-forward distributions H of Ji, so we start 
with the evolution equations [[| for H(x,£) with 16 [—1, 1] and use transformations fl2T|) 
and (|24]) to rewrite them in terms of the distributions F{X, () with X G [0, 1]. 

In the DGLAP-like region X > ( (which corresponds to x > £ or x < — £) the equations 
that we obtain for JF are equivalent to those given for the non-forward distributions of 
Radyushkin [EL M. Their full form can be found in the Appendix. Moreover in the limit 
C — > they reduce to the familiar DGLAP evolution equations. 

However in the ERBL-like region X < ( (corresponding to — £ < x < £) the equations 
obtained for J 7 are different to those given in || |9| for the non- forward distributions. 
They have the following forms 



+ 



Pc 



v 



a s C F r 1 dZ 
k ~Z~ 



71 



x-c + z 



X 



^(Z,() 



(29) 



+ 



ocsCp f 1 dZ 

k ~z~ 



Pi 



7T 



To 



asNf 



7T 



Z 



X 

T 



x-C + z 

idZ (1-C/2)(C-X) 



c z 



f s (z,() 



AX 2X - 

T + — 



(30) 



+ 



Pgq ® F b + Pgg ® 



a 5 C F /-idZ (C-Xf 



7T 



c ^ C(l - C/2) 
1 dZ (C - A) 2 



7T 



1 



2Z / 2A X 



(31) 



where the scale [i is implicit in the distributions f . The full forms of the equations 
are given in the Appendix. Here it is sufficient to note that the convolutions shown 
symbolically as P <g) T are identical to those given in [pj, |9| . However the new evolution 
equations contain several additional terms, each being a convolution integral over the 
range [(, 1]. These extra terms are essential to preserve the symmetry properties ( p?D and 
(|2Sj) of T during the evolution. We note that in the limit £ — ► 1 the additional terms are 
to equal zero and that ( |2"9"D - fl3~T| ) reduce to the ERBL evolution equations |T2[ for the 
distribution amplitudes. 



8 



5.1 Numerical results of the evolution 



To illustrate how the off-diagonal distributions T evolve with increasing renormalization 
scale /i we constructed a computer programme based on the equations given in the Ap- 
pendix. For the initial input at the starting scale \i = 1 GeV we adopt the following 
strategy. We start with given input forms for the off-forward distributions H gtg (x,£), 
which are even in £. An example for the quark distribution is shown in Fig. 5(a). Then 
using prescriptions (|2l|) and (p4j) we transform H q>g (x,^) into the distributions Fq,q, g (X, () 
which satisfy the symmetry relations ( ^7|) and (^8|). The initial distribution H q (x,!;) shown 
in Fig. 5(a) is only meant to illustrate the general features of the adopted strategy. The 
detailed properties of more realistic initial distributions will be discussed in a separate 
paper. 

The results that are obtained by evolving jF y , JF 5 and T g to higher scales are shown 
in the three plots of Fig. 6. In each plot the dashed curve is the input at fi = 1 GeV, while 
the dot-dashed curve shows the effect of evolution up to n = 10 GeV. It is evident that 
evolution does indeed preserve the symmetry properties in the ERBL-like region, X < 

The continuous curves in Fig. 6 are the results of evolving all the way up to — > oo. 
These asymptotic forms are identical with the analytic asymptotic solutions of the 
evolution equations for the distributions T given in the Appendix 

f s (X,() ~ -(l--)l — -I) (32) 




WO ~ 

A remarkable property || is evident from Fig. 6. We see that the distributions are swept 
from the DGLAP-like to the ERBL-like region as fi increases. Indeed the asymptotic 
forms show that they are finally entirely contained in the ERBL-like region with X < 



6 Relation to the non-forward distributions 

The off-diagonal distributions F(X, (), constructed in the previous section, are equivalent 
to the off-forward distributions H(x,£) defined by Ji. They are also closely related to, 
but not the same as, the non- forward distributions J~c,{X) introduced by RadyushkinF]. 
The difference between them occurs in the ERBL-like region (X < (). 

5 We thank A.V. Radyushkin for helpful comments on the subject of this section. 
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The non-forward distributions ^^(X) are related to the off- forward distributions 
H q (x,£) in the following way (see Section IX of Ref. || for a detailed discussion) 



[ ^(X) if x > £ 

^(X)-^(C-X) if-£<x<£ (33) 

-?!{{- x) ifx<-e, 



where X = (x+0/(l+0 and C = 2£/(l+£). Notice that while in the DGLAP-like regions 
(x > £ or x < — £) there is a one-to-one correspondence between the two distributions, 
in the ERBL-like region (— £ < a; < £) Ji's distribution H q only determines a specific 
combination of Radyushkin's distributions T q,q . This is in contrast to the distributions 
defined by eqs. (|21~D which are in one-to-one correspondence with H q . 

Comparing eqs. ( |33"D with eqs. (|21"D we see that the off-diagonal distributions T are 
identical to the non-forward distributions T in the DGLAP-like region [X > £). However 
in the ERBL-like region there are different. To be precise, for X < (, we have 

^(x,c) = ^(x) - ^((-x) 

(34) 

^(x,c) = J* c {x) - ^(C-x). 

The main difference between the distributions !F{X, Q and the non-forward distributions 
JF^(X) is that the latter do not obey the symmetry properties (|25| ) and (p?|)-(p^). These 
properties are essential for our distributions. They result from the construction which 
ensures the equivalence of the distributions jF to Ji's distributions H. The physical 
reason for the symmetries was discussed in Section 4.2. An important consequence of the 
symmetry relations is that in the ERBL-like region the quark and antiquark off-diagonal 
distributions are not independent, see relation (f25|). 

This should be contrasted to the case of the non-forward distributions of Radyushkin. 
They are are obtained through the integration of "double distributions" F which are 
universal (^-independent functions. The double distributions are separated into two inde- 
pendent components (which are denoted by F q and F q ) according to the sign of x in the 
exponential. As a result the corresponding non-forward distributions T q and T q are also 
independent in the ERBL-like region, see for more details. Thus there are twice as 
many quark "degrees of freedom" in the ERBL-like region as in our case. 

A similar comparison can be done for the non-singlet, singlet and gluon distributions. 
As a result we find the following relations for X < £ 

T V {XX) = J¥(X)+J%({-X) 

P S (X,C) = J-f(X)- J-f(C-X) (35) 

P(x,c) = ^(x) + ^(c-x). 
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Thus we see that the distributions T are equal to symmetric or antisymmetric combina- 
tions of the corresponding non-forward distributions T in the ERBL-like region. These 
combinations of the non-forward distributions were used in || in the description of the 
ERBL-like region. However the further analysis in Ref. || was done in terms of the 
unsymmetrized non-forward distributions J-j- V ' S,9 \X). 



6.1 Comparison of the two sets of evolution equations 

The evolution equations for the non-forward distributions JF^' 5,51 ) of Refs. ]5|, |S| do not 
obey the symmetry properties (|27j) and (^) in the ERBL- like region. One may try, 
however, to write down the evolution equations for the combinations on the right hand 
side of eqs. (p5|), starting from the equations given in H |9] for the full non-forward 
distributions 

jr^x) EE ^ sym) (X) + jf sym) (X) 

(36) 

= | [F C {X) + ^ c (c -x)\ + | [J- C (X) - ^ c (c - X)} 

Not surprisingly these "symmetrized" evolution equations are almost identical to the 
evolution equations for the distributions (p9|)-([3l|) . The integrals over [C,l], indicated 
explicitly in fl29|)-(|31|), appear in the "symmetrized" equations as a result of the sym- 
metrization procedure. The only difference appears in the "symmetrized" gluon equation 
which additionally contains a term proportional to the integral over the full non-forward 
singlet distribution 

£dZ ?S(Z) = £ dZ jf sym) (Z) + J^dZ^(Z) . (37) 

Since the "symmetrized" gluon equation should contain only the asymmetric singlet com- 
bination ^ asym ^ ; the above term mixes the symmetric and antisymmetric components of 
the singlet distribution, and thus violates the symmetry properties (|27 ) and (p8) for the 



singlet and gluon distributions. The only case when it does not happen is if the integral 
( 371) is equal to zero due to initial conditions. The value of this integral is conserved by 



the evolution equations [|5], [| for the full non-forward distributions^. Only in this case 
may the non-forward distributions of Radyushkin be equivalent to the off-diagonal distri- 
butions of Ji. This can be done by taking into account only one of the two parts in the 



decomposition ( p6|) — symmetric for the non- singlet and gluon, and antisymmetric for 
the singlet, distributions. 



6 We were informed by A.V. Radyushkin that the above mentioned problem with the integral (g7|) can 
be solved if one uses the kernel Pgq in the evolution equations of Ref. || in the form originally obtained 
by Chase @. The method used in Ref. [| cannot unambiguously fix this kernel. 
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6.2 The singularity structure of the basic amplitude 

For the purpose of illustration we may consider the classic process of deeply virtual Comp- 
ton scattering. The invariant amplitude for the process has the generic form |2| 

T ~ J dx 

If the amplitude is translated into a form involving the distributions T q (X, Q defined on 
the interval < X < 1, then fl38|) becomes 

T „ / dX + + f « [Ux , () + MXi 0] . (39) 

Jo X — Q + is Jc, X 1 J 

We see that ( |3"9"D contains only one singularity at X = (, which results from the quark 
propagator, and is regularized by the +ie prescription and assuming that F q>q (X, () are 
continuous at X — £. Note that there is no singularity at X = since the second integral 
is bounded by ( > from below. 

This is in contrast to the amplitude derived using the non- forward distributions T^{X) 
0, H]. Then T contains a second singularity at X = 0, since in this case the second 
integral in (|39|) goes from to 1. The result can be derived by substituting relations 
fl34] ) into ([39]). This additional (end-point) singularity is removed by assuming that the 
non-forward distributions T^' q {X) vanish as X —>■ 0. Looking at eq. (|33|) we see that this 
assumption is equivalent to the continuity of H(x,£) at x = ±£ (or JF(X, £) at X = 
Such assumption is not required for our off-diagonal distributions jF(X,Q, see eq. fl3"9"D- 
Indeed, if present, it would clearly violate their continuity at X = (, or their symmetry 
about X = (/2, see Fig. 5. 



+ 



x 



£ + ie x + £ — is 



(38) 



7 Conclusions 

In this paper we have transformed the off-forward parton distributions H(x,£) defined 
by Ji, in which the defining direction is the average between the incoming and outgoing 
proton momenta and x G [—1,1], into off-diagonal distributions F(X, (), in which the 
defining direction is the incoming proton momentum and X G [0, 1]. These off-diagonal 
distributions f{X, Q therefore have a close identification with conventional (diagonal) 
distributions. Moreover, by construction, they are fully equivalent to the off-forward 
distributions of Ji. 

In the ERBL-like domain (X < () they satisfy the symmetry relations 

HC-xx) = ±H*,Q (40) 

where the + sign applies to the gluon and quark non-singlet distributions, and the — sign 
applies to the quark singlet. We presented the evolution equations satisfied by the JF{X, () 

7 We thank A.V. Radyushkin for this remark 
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and gave numerical results (Figs. 5 and 6) to illustrate the properties of the distributions. 
We found that asymptotically (// — > oo) the distributions evolve to the known analytic 
asymptotic forms. Indeed as \x increases the distributions are swept from the DGLAP-like 
domain to lie entirely within the ERBL-like region, as illustrated by the example shown 
in Fig. 6. The symmetry relations ([It]) are preserved at each stage of the evolution. 

The distributions T{X, Q are analogous to, but not the same as, the non-forward 
distributions JF^(X) introduced by Radyushkin ||]. The difference lies in the ERBL-like 
region, since the non-forward distributions do not obey the symmetry relations (|40|). As a 
result the non- forward distributions JF^(X) are not in general equivalent to the off- forward 
distributions H(x,£) of Ji. We stressed that this happens only in the ERBL-like region. 
We discussed conditions under which T^{X) would become equivalent to H(x,£) (and 
F(X, 0). We also commented on the singularity at X = of the basic DVCS amplitude 
at tree level when written in terms of JF^(X), which requires T^{X) to vanish as X — > 0. 
The distributions F(X, (), which we defined, have the advantage that they do not lead 
to such a singularity. 
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Appendix 



Here we present for reference the full form of the evolution equations for our non-singlet 
!FY(X, C, /j,), singlet f 7 (X, (,fi) and gluon T g {X, distributions defined in the range 
< X < 1 by eqs. ( p6[ ) and (|24]). The asymmetry parameter ( lies in the range [0, 1]. 

We use the following notation X' = X — ( and Z' = Z — Q and suppress the renor- 
malization scale /i among the arguments of our distributions. In the DGLAP-like region 
X > ( we have the following evolution equations 








^r s {xx^) 



a s 



71 



c h 



dZ 



■ x 
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X X' 



+ f s (x,() 



z z 
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dZ 
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Yz 1 
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[(x/z) + (x'/z')] ^(x,c) - U/zf + (x'/z') 2 ] ^(z,c) 



+ -WO 



X-Z 

11 - (2ty)/3 + ln (1-X) 2 



2X. 



i-C 



(41) 



where Cp = 4/3 and N c = 3, and X/ is the number of active flavours. In the limit ( = 
the above equations become the familiar DGLAP evolution equations. 

The equations in the ERBL-like region X < ( are more complicated since they involve 
integration with different kernels in the intervals [0,X] and [X, 1]. We have 
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X' 



t v Azx) , fY{z,0-f s {x,Q- 



C 



+ 



x-z 



+ 



W0 + 5(^)-^ft() 



+ A y (x,c,/i) 



C 

3 , X(l-X) 



Z-X 



+ 



1 dZ_ 
c 2~ 



+ 







c 



+ 



i 



/ 



j*(z,c) , ^(z,c)-^(x,c) 



+ 



1 dZ 



3 , 
2 +ln 



C + z-x 

X(l-X)" 



C 



X 





N f 


7T 




1 


dZ 


/ 










+ 


r 



x-(+z c 
dz, c\ ,x 

C 2 V 2) \Z> 



X 2X - C 
4— + 



C C-z 



X\ C-2X 

in -c) + V- 



4X 2X - C' 



WO 



d - 



a s 



IT 



dZ 



'z 7 



1 - 



CJ i-C/2 



+ 



jdz(. 



x*\ P^zx) 

Z(J l-C/2 



idz (c-x) ; 



Jc ~~Z 



C(i-C/2) 



f s (z,() 



+ ^x c 

7T 



/ 



x^ 



X 



X 

1 + 2— + 



X 



C C-2 



r 9 (z,() 



16 



+ 



(x'/z')f g (z,()-f 9 (x,0 
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(42) 



+ 



i dZ (c - x) s 
c ~~z z 



1 



2Z ( 2X X 

x-c + z + cn 1 + T + ^, 



^0 



For C — 1 the above equations reduce to the ERBL evolution equations for the distribution 
amplitudes. It is also instructive to check that both set of equations, ([41]) and (f42|), lead 
to the same limiting set of equations when X — > ( from both sides. 

The equations for the singlet JF S and the gluon T 9 distributions form a coupled set 
of equations which, in general, need to be solved simultaneously in both the ERBL- and 
DGLAP-like regions. However for X > ( it is sufficient to solve the equations only in 
the DGLAP-like region since the integration in (IT) involves only parton distributions for 
values of Z > X (as is true for the DGLAP equations in the limit ( = 0). This is not 
the case if X < £. Then the solutions depend on the values of the parton distributions in 
the full interval [0, 1], and so both the set of equations, ([41]) and (0), have to be solved 
simultaneously. 
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(a) x>0: q(x) = H q (x) 



(b) x<0: q(-x) = -H q (x) 



b V 



d d 





Figure 1: Schematic diagrams showing the contributions to H q (x) with x > and x < 
respectively, which can be identified with the familiar quark and antiquark distributions. 
b, are the quark annihilation and creation operators and d, d^ are those for the antiquark. 
The momentum fractions refer to the plus light-cone component of the incoming proton 
momentum P. 
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(a) x>h;. DGLAP-type region for 
the quark distribution 




(b) -^<x<^: ERBL-type probability 
amplitude 




(c) x<-L,\ DGLAP-type region for 
the antiquark distribution 




Figure 2: Schematic diagrams of the off-diagonal distribution H q {x,£), in the three 
distinct kinematic regions. The proton and quark momentum fractions refer to P + , where 
P is the average of the incoming and outgoing proton four momentum. Note that the 
four momentum transfer satisfies A + = 2£P + and that x covers the interval [—1,1]. 



19 



Figure 3: The proton and quark momentum fractions with respect to the initial pro- 
ton momentum P corresponding to the off-diagonal distributions F(X, () defined in the 
domain < X < 1. The four momentum transfer satisfies A + = (P + . 




Figure 4: A sketch showing how the support — 1 < x < 1 of the off-diagonal distribution 
H q is translated into the regions < X < 1 of the two functions T q and Tq. The 
translations are given by (|T^) and (pip, or by the inverse relations fl2"2"|). 
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Figure 5: (a) An example of the off-diagonal distribution H q (x,£) with £ = 0.5; (b) the 
distributions F q (X, () and Tq{X, Q generated from H q (x,£), and (c) the resulting non- 
singlet jFq and singlet jF s distributions showing their symmetry and antisymmetry in the 
ERBL-like region X < (. 
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Figure 6: Evolution of the non-singlet , singlet T s and gluon T g distributions defined 
in the range [0,1] from initial input at \i — 1 GeV (dashed curves). The asymmetry 
parameter ( = 0.5. The dotted and continuous curves correspond to fi = 10 GeV and 
fi — > oo respectively. The latter curves are identical to the analytic asymptotic solutions 
given in (|32j). 
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